APS/123-QED 



OO 
O 

o 

(N 
O 

Q 

(N 



CZ3 



Heterogeneous nucleation and adatom detachment at 1-D growth of Indium on 

Si(100)-2xl 

Jakub Javorsky, Martin Setvm, Ivan Ost'adalQ and Pavel Sobotfk 
Department of Surface and Plasma Science, Charles University, 
V Holesovickdch 2, 180 00 Prague 8, Czech Republic. 

Miroslav Kotrla 
Institute of Physics, Academy of Sciences of the Czech Republic, 
Na Slovance 2, 182 21 Prague 8, Czech Republic 
(Dated: December 29, 2008) 

Growth of atomic indium chains - ID islands - on the Si(100)-2xl surface was observed by 
scanning tunneling microscopy (STM) at room temperature and simulated by means of a kinetic 
Monte Carlo method. Density of indium islands and island size distribution were obtained for various 
deposition rates and coverage. STM observation of growth during deposition of indium provided 
information on growth kinetics and relaxation of grown layers. Important role of C-type defects at 
adsorption of metal atoms was observed. Measured growth characteristics were simulated using a 
microscopic model with anisotropic surface diffusion and forbidden zones along the metal chains. 
An analysis of experimental and simulation data shows that detachment of indium adatoms from 
the chains substantially influences a growth scenario and results in monotonously decreasing chain 
length distribution function at low coverage. Diffusion barriers determined from the simulations 
correspond to almost isotropic diffusion of indium adatoms on the surface. The results are discussed 
with respect to data reported in earlier papers for other metals. 

PACS numbers: 68.55.A-, 68.55.ag, 68.37.Ef, 81.15.Aa 
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I. INTRODUCTION 



The Si(100)-2xl surface is composed of silicon atom 
pairs - dimers - arranged into rows. It represents a 
natural template for spontaneous growing linear struc- 
tures of many materials, like group III-V metals. Group 
III metals (Al, Ga, In, Tl) are known to grow in onc- 
dimensional (ID) atomic chains when deposited on the 
Si(100) surface 1 ^^. Technique of scanning tunneling 
microscopy (STM) enabled detailed study of metal lay- 
ers with atomic resolution. An especially powerful tool 
for investigation of growth kinetics is the STM in-vivo 
technique 5 ^ which allows direct monitoring of the layer 
as it grows. ID metal chains grow perpendicularly to the 
underlying silicon dimer rows of the Si(100)-2xl surface. 
They are composed of metal dimers (oriented parallel to 
the silicon dimers) sitting in the "trenches" on the silicon 
surface 7 . Growth of the metal chains has been explained 
by a surface polymerization reaction^. Metal chain ends 
act as nucleation centers. Since the sites adjacent to a 
chain are energetically unfavorable for adsorption (no ad- 
sorption has been observed there) the chains grow only 
in length. The chains are always separated by a dis- 
tance of at least 2a (a= 0.384 nm, surface unit cell spac- 
ing). Thus, the surface is saturated by metal adatoms 
at a coverage of 0.5 ML (1 ML= 6.78 x 10 14 cm" 2 ). 
Some differences exist between various group III metals. 
While Al and Ga chains are believed to be stable at room 
temperature^, indium atoms are known to detach from 
chains and re-attach to other chains^ and Tl chains were 
shown to be even more unstable^. Though a qualitative 
description of diffusion and growth processes for group 



III metals exists, the corresponding values of microscopic 
parameters are not known. The heights of diffusion bar- 
riers on Si(100) have not been yet determined for In and 
the values reported by Albao et al. 9 for Ga and by Brocks 
et al& for Al are very different. Similar discrepancy ex- 
ists for estimation of Ga dimer pair-interaction energy - 
Tokar and Dreysseii suggest « 0.2 eV while Takeuchi's 
ab-initio calculation gave rj 0.8 eVi2. Recently Kocan 
et al. reported 1 ^ that detachment of indium atoms from 
chains is length dependent, so interactions other than 
nearest-neighbor (NN) might play a role. 

The role of surface defects present on the Si(100)-2xl 
surface 1 ^ at the metal adsorption and nucleation was 
reported and discusse d 10 : 14 ' 15 : 16 . Experimental results 
showed that influence of A- and B-type defects (one and 
two missing dimers, respectively) can be neglected but 
C-type defects are important. The C-type defects, which 
appear on STM images as a small bright protrusion next 
to a dark spot in filled states and as a larger bright spot 
in empty states, are reported to be very reactive and act 
as nucleation center a 10 ' 14 . The C-defects were indepen- 
dently interpreted by Hossain et al*^ and Okano and 
Oshiyamai^ as dissociated H2O molecules, with the H 
and hydroxyl group bonded to neighboring silicon atoms 
of two adjacent surface dimers. The results both of ex- 
perimental and theoretical study of In nucleation at the 
C-defects were reported in Ref. [l6|. After adsorption of 
an In adatom at a C-defect (exclusively on the unoccu- 
pied side of the two adjacent Si dimers) the chain begins 
to grow in one direction only (see Fig. 1 in Ref. 16). The 
chain termination at the defect is stable, the opposite 
"free" end is active as an adsorption site for adatoms. 



Albao et al£ reported growth characteristics for Ga 
on Si(100) at room temperature (RT). An unconventional 
monotonously decreasing scaled island (chain) size distri- 
bution function obtained for low coverage was explained 
by an irreversible growth model and kinetic Monte Carlo 
(KMC) simulations. The simulations resulted in the 
monotonous size distribution only if highly anisotropic 
surface diffusion of Ga adatoms was introduced other- 
wise a monomodal form of the distribution function was 
obtained. The presence of C-defects included late r 14 ! 15 
did not change the results significantly. Similar growth 
characteristics at RT we reported for indium 1 ^. The 
chain length distributions obtained for various coverages 
are monotonously decreasing and obey a scaling relation. 
Most of In chains (60-90 %) in observed layers were on 
at least one end terminated by a C-type defect. In the 
in-vivo experiments the percentage was higher, 90-100 
% (due to low deposition rates of In during the in-vivo 
measurements). Density of indium chains and average 
chain length depend on C-defect concentration. Another 
phenomenon (observed at higher coverage) is that indium 
atoms are able to migrate throughout sites adjacent to 
an indium chain even though no adsorption is observed 
(such an effect seems to be negligible at low coverage). 

In this paper, we use STM data and KMC simulations 
for detail studying indium growth on the reconstructed 
surface Si(100)-2xl. A growth model, which includes 
presence of C-defects and a process of atom detachment 
from indium chains, is used for studying a role of C- 
defects at metal chain nucleation, determination of diffu- 
sion barriers and investigation of relation between growth 
processes and a form of chain length distribution func- 
tion. Processes and parameters included in the growth 
model are discussed with respect to experimental data 
obtained by means of STM. 

Performed experiments are characterized in Sec. 
II, consequently experimental results are presented in 
Sec. III. A simulation model is described in Sec. IV, re- 
sults of simulations are compared with the experimen- 
tal data and discussed in Sec. V, formulas for calcula- 
tion of deviations between experimental and simulated 
growth characteristics are given in the Appendix. Fi- 
nally, Sec. VI contains our conclusions. 



II. EXPERIMENTAL 

All experiments were carried out at room temperature 
in a non-commercial UHV STM system with base pres- 
sure < 3 x 10~ 9 Pa. Si(100) samples were cut from an 
n-type, Sb doped silicon wafer with resistivity < 0.014 
57cm. To obtain the 2x1 reconstruction, samples were 
flashed several times for « 20 s to « 1200° C. Indium 
was deposited from tungsten wire evaporators either be- 
fore or during imaging the surface by means of STM. In 
the latter case (experiment in-vivo) a miniature evapora- 
tor was in a distance of 4 cm from the sample and a beam 
of In atoms was determined by means of two apertures 
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FIG. 1: Filled state STM image of indium ID chains grown 
at room temperature on the surface Si(100)-2xl. The bright 
spots at ends of some chains indicate termination by a single 
indium ato m 10 ' 16 ; stable terminations at C-defects are marked 
by the arrows perpendicular to chains. Unoccupied C-defects 
are marked by the arrows parallel with chains - the arrows 
are directed from the side of possible chain growing. Cover- 
age 0.044 ML, in-situ measurements 2 hours after deposition; 
image size 30x30 nm 2 . 



with a diameter of 1 mm. Incidental angle of the beam 
was 30°. The apex shape of a tungsten electrochemically 
etched tip enabled deposition "under" the STM tip onto 
the scanned area. The thermal drift during deposition 
was compensated by the STM control unit. The in-vivo 
measurements provided continuous STM imaging of the 
investigated surface area before and during the deposi- 
tion at a rate of 1 image/min. At standard in-situ mea- 
surements the indium layers were observed 0.5-4 hours 
after deposition. We used a tip voltage of +2 V and tun- 
neling current 0.3 nA, the values at which tip influence 
on the detachment of indium adatoms from the chains is 
minimized as we proved earlier—. 



III. EXPERIMENTAL RESULTS 

An example of indium chains grown on the Si(100) 
surface at low coverage is shown in Fig. [T] Length of 
the chains can be easily determined with atomic preci- 
sion from filled states images where the "free" termina- 
tions of chains containing an odd number of atoms appear 
much brighter than in case of even numbe r 10 ' 14 . Chain 
terminations at C-defects and unoccupied C-defect are 
marked by arrows. Concentration of the C-defects in- 
creases during deposition at in-vivo experiments (proba- 
bly due to thermal desorption of residual water molecules 
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FIG. 2: Dependence of average chain length on coverage ob- 
tained from STM experiments. The average length of In 
chains increases with coverage and saturates at a coverage 
of 0.25 ML for the in-vivo data. The dependence increases 
faster for the layers prepared at in-situ experiments due to 
smaller and constant concentration of C-defects. The effec- 
tive deposition rate for C-defects was 4xl0 -6 ML/s during 
the in-vivo measurements. 



from heated parts near the evaporation source in a rel- 
atively small distance from the sample). The increase 
was found linear and corresponds to a deposition rate 
of 4 x 10~ 6 ML/s. The deposition rate of C-defects at 
in-situ experiments was negligible (due to a large dis- 
tance of the evaporation source and better screening of 
the sample) and the initial concentration of the C-defects 
can be considered as unchanged. 

Since the atomic structure of In chains is well 
know n 7 ' 10 ' 16 we focused on statistical characteristics of 
the chains. Obtained data are summarized in the Ta- 
ble HI The data were acquired from images of the size 
of 30x30 nm 2 or 40x40 nm 2 . The image areas are large 
enough (30 nm corresponds to a ~ 100-atoms long chain) 
and the resolution is sufficient to discern the number of 
atoms in the chain. Data were collected only from ter- 
races much wider than an average chain length to ex- 
clude the influence of step edges on the chain growth. 
Statistical characteristics were evaluated for both types 
of growth experiments (in-vivo and in-situ). 

Length of the grown indium chains evolves in time due 
to attachment/detachment of single atoms to/from the 
"free" ends. Fig. [2] contains dependence of an aver- 
age chain length on coverage obtained from two in-vivo 
measurements with deposition rates 6xl0 -5 and lxlO -4 
ML/s and a number of in-situ measurements for var- 
ious deposition rates (from 2xl0 -3 to 3xl0~ 2 ML/s) 
and coverages (0.01-0.15 ML). The average length of 
chains is smaller in case of in-vivo experiments because of 
higher concentration of C-defects, and from the coverage 
0.25 ML (which corresponds to occupation of a half of all 
possible adsorption sites on the surface) almost does not 
increase. 



FIG. 3: Histograms of chain length distribution obtained 
for two indium layers with coverage of 0.04 and 0.08 ML 
deposited at a deposition rate of «0.004 and «0.002 ML/s 
respectively (in-situ experiments). For the comparison his- 
tograms contain single indium atoms trapped on C-defects. 
N represents a number of investigated chains. 

Fig. [3] shows histograms of chain length distributions 
for two different coverages obtained from in-situ experi- 
ments. The data were averaged over STM images taken 
after the deposition from various surface areas. The 
chains nucleated on C-defects are distinguished from the 
"free" chains (without termination on a C-defect). The 
histograms contain also single In atoms adsorbed and 
trapped on C-defects. They represent stable objects with 
a role of nucleation centers. Their presence in the his- 
tograms allows better understanding of the growth but 
we consider only chains with length s > 2 as parts of 
the investigated "island population" . The histograms 
are decreasing for s > 1 and the same tendency was ob- 
served for the other prepared layers with low coverages 
< 0.15 ML. The monotonously decreasing chain length 
distributions obtained for growth of indium at RT rep- 
resent a remarkable quality. It can be seen that his- 
tograms contain some features related to details of the 
chain growth - a small excess of chains containing even 
number of atoms which corresponds to higher stability of 
such chains experimentally observed^ - and calculated^* 
before. 

The chain length distributions exhibit scalin g 20 ' 21 : the 
function N s (s) 2 /Q scales with s/(s), where N s is den- 
sity per site of chains of the length s, (s) - average chain 
length and 8 - coverage. The upper panel (a) of Fig. [4] 
shows the scaled distribution functions corresponding to 
various delay between the end of deposition and STM 
measurements at the in-situ experiments. Only chains 
with s > 2 were included. The all data were fitted by 
an exponential function. Due to statistical fluctuations 
in the data it is difficult to distinguish reliably an effect 
of postdeposition relaxation. The postdeposition relax- 
ation can be expected because of the process of detach- 
ment, which introduces a feature of "reversibility" into 
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TABLE I: Measured statistical characteristics of In chains for various coverages and deposition rates. 
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a growth mechanism. We will discuss the "reversibility" 
and postdeposition relaxation later in Section V. 

The scaled distributions obtained from the images 
recorded at in-vivo measurements are in the bottom di- 
agram (b) in Fig. [4] The data suffer from limited size of 
the investigated surface area and relatively small num- 
ber of observed metal chains. Values for a particular dis- 
tribution corresponding to a chosen moment (coverage) 
of the growth were averaged from a set of 3 subsequent 
STM images around that moment. The whole recorded 
sequence of images taken during the growth up to 0.2 ML 
at a deposition rate of 0.0001 ML/s (estimated from the 
images) is represented " equidistantly" by the 15 distri- 
butions. The rather scattered data exhibit monotonous 
character and can be approximated by an exponential 
function. 



IV. SIMULATION MODEL 

We developed a physical atomistic model for submono- 
layer growth of indium on Si(100). Similarly as in the 
model used by Albao et al. to describe growth of gal- 
lium layers at low coverage^, we consider anisotropic dif- 
fusion of the metal adatoms on the 2x1 reconstructed 
surface. However, we took into account experimental re- 
sults obtained for In and included new features: presence 
of C-defects acting as nucleation centers on the surface 
and detachment of single atoms from chains. The mech- 
anism of detachment introduces into the growth model 
a possibility of "reversible" behavior ("reversible" does 
not mean here "equilibrium" - due to existing flux of 
deposited atoms). Approaching an equilibrium state de- 
pends on growth conditions and it competes with kinetics 
which controls growth entirely in case of an irreversible 
model used for gallium. 

The model for In deals with four types of objects on 
the surface (see Fig. [5]): 

1) Free In adatoms - indium atoms perform thermally 
activated hopping on the surface represented by a square 
lattice (see Fig. [5]). The hopping is anisotropic (par- 
allel and perpendicular to indium chains) described by 
rates h n = vexp[— E\\ ± /(fesT 1 )], where the attempt fre- 
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FIG. 4: Scaled island size distribution functions obtained 
from experimental data, (a) Results from three in-situ mea- 
surements performed at various moments after deposition, (b) 
Data obtained during deposition (in-vivo) - composed from 
a set of 15 histograms covering equally the whole deposition. 
The data were fitted by exponential functions. 



quency was set as v = 10 13 s _1 , the activation energies 
(diffusion barriers) E\\ , E± are simulation parameters, 
ks~ the Boltzmann constant and T- temperature. 

2) Forbidden zones - surface sites not accessible for 
hopping adatoms. The forbidden zones were introduced 
- similarly as in Ref. 9 - to simulate the ID growth, min- 
imum separation observed between two adjacent chains 
and the fact that no chain nucleates at a "non-reactive" 
side of a C-defect. 
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FIG. 5: Left: Filled state image of indium chains. The two 
neighboring chains are separated by a distance of one lat- 
tice unit. Dark spots at the ends of chains are C-defects. 
Right: Alignment of the model to Si(100) substrate. The 
black square represents the "reactive site" at the C-defect 
position, grey squares - adsorption sites occupied by in- 
dium adatoms, empty squares - unoccupied adsorption sites, 
crossed squares - "forbidden zones"- energetically unfavor- 
able adsorption sites. 



3) C-type defects - dissociated H2O molecules - rest 
on top of dimer rows. STM observations show that one 
side of the C-defect acts as a nucleation center while no 
adatom adsorption has been observed on the opposite 
"dark" side (on filled state STM image) of the defect. 
Thus, C-type defects are represented by a single forbid- 
den site and an adjacent "reactive" site. 

4) Indium chains - orientated perpendicularly to the 
Si dimer rows - are composed of two or more In atoms. 
Similarly as in the model used by Albao et al£ arrange- 
ment of atoms in chains into dimcrs is not taken into 
account. 

There are three main processes in our model: 

A) Deposition - Adatoms and C-defects are deposited 
randomly. If a defect or an adatom is deposited on an 
already occupied position or (in case of In atoms) into 
a forbidden zone, the nearest free position is looked up 
in the direction of dimer rows. Deposition flux and time 
were set the same as in the particular experiments (in- 
dium flux Kr 2 -lCr 5 ML/s, C-defect flux 4 xlCT 6 ML/s 
- for "in vivo" only) and corresponding simulations were 
performed for all STM experiments. Orientation of a de- 
posited C-defect (reactive site) was chosen randomly. As 
the C-defects change their state only very rarel y 16 ' 22 the 
orientation is fixed during the whole simulation (it differs 
from Ref. fill where the orientation was determined each 
time when an adatom tried to hop next to the defect). 

B) Surface migration (hopping) - Single indium atoms 
deposited on surface perform random hops among ad- 
sorption sites. Jumps to forbidden zones and on top of 
other atoms are prohibited An atom is trapped when 
hops onto a C-defect reactive site. If an atom hops on a 
site next to another atom in the direction parallel to In 
chains (perpendicular to the dimer rows), a new indium 
chain nucleates or an existing one grows. The hopping 
rates h\\ t ± in directions parallel or perpendicular to the 



indium chains are given by the simulation parameters E11 
and Ex (diffusion barriers parallel and perpendicular to 
indium chains). 

C) Detachment - An indium atom can detach from a 
chain end, not terminated by a C-defect, by jumping off 
either in perpendicular or parallel direction to the chain. 
According to our best knowledge there are neither exper- 
imental nor theoretical data available to characterize the 
detachment direction in the studied (or similar) system 
on the atomic level. The process of the detachment is 
thermally activated and can be described by two param- 
eters: -E(Zei|| and Ed e t± (energy barriers for detachment 
in the parallel and perpendicular direction). Our previ- 
ous measurements^ , and theoretical calculations^ 6 - show 
that the energy for detachment depends on length and 
termination of an In chain by a single atom or dimer. 
The detachment from a chain containing an odd number 
of atoms is easier than in case of even number. For the 
simplicity the model contains only one parameter - the 
value of a mean energy barrier for detachment, 0.82 eV, 
derived from experimental data reported in Ref. [l(| As 
we know only the total rate of detachment without any 
details, simulations were performed for detachment cither 
in the parallel or perpendicular direction. Two-atomic 
chains (dimers) on C-defects are considered as stable ob- 
jects which cannot decay (both experiment and theoret- 
ical calculations^ 6 - confirmed their high stability). 

The simulation proceeds as follows: C-defects are ran- 
domly deposited on the surface with initial concentration. 
Then indium atoms are randomly deposited (together 
with additional C-defects according to a simulated exper- 
iment). For details of an employed method of KMC sim- 
ulations see Ref. [23|. After the deposition stops the layer 
is allowed to relax for the same time as in the correspond- 
ing experiment. The Si(100)-2xl surface is represented 
by a square matrix of adsorption sites, each capable of 
holding a single indium atom. We used a matrix size 
between 100x100 and 500x500 lattice units, each sim- 
ulation was repeated at least 9 times and the obtained 
data were averaged. The values of the matrix size were 
chosen so as the mean average error of simulated data 
was 2-3 x smaller than that of measured data (size of 
statistical arrays varied for different values of coverage) . 
The boundary conditions were periodic. 



V. RESULTS OF SIMULATIONS AND 
DISCUSSION 

Diffusion barriers. The in-situ and in-vivo experi- 
ments were simulated for various combinations of diffu- 
sion barriers Eh , E± . Comparison of experimental and 
simulated growth characteristics (average chain length, 
dependence of average chain length on coverage and 
scaled chain length distribution) was used for estima- 
tion of a combination of barriers which provides the best 
agreement - the lowest deviation calculated as presented 
in the Appendix. 
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a) Average chain length b) Average length of chains c) Chain length d) Average chain length 

in-situ not terminated by C-defects distribution function in-vivo 




£[eV] £„[eV] E a [eV] £ [eV] 

FIG. 6: Maps showing deviation between the measured and simulated characteristics of indium chains as a function of diffusion 
barriers parallel and perpendicular to In chains. White areas indicate the best agreement, black areas the worst. Dotted line 
connects local minima in panel a), dashed line in panel b). Dashed line from panel b) is also shown in panel a) for comparison. 
Energy barrier for atom detachment was 0.82 eV. Only detachment parallel to indium chains was allowed in these simulations. 



Figure [5] contains diagrams with the dependence of 
the calculated deviations on the simulation parameters 
E\\ and E± for the chosen growth characteristics. The 
grayscale represents accuracy of the match for a given 
combination of parameters with white for the most pre- 
cise match. The plots a) and b) demonstrate fitting the 
in-situ measurement for a layer with coverage 0.08 ML. 
In the plot a) the dotted black line shows the combina- 
tions of activation energies which provide the best match 
between the measured and simulated average length of 
indium chains, the dashed line in the plot b) shows the 
best agreement for the average length only of those chains 
which arc not terminated by C-defects at any of ends. 
The relatively small "sub-population" of these chains be- 
haves in a different way so the two plots together were 
used to determine the best combination of energies: Eu = 
0.62±0.03 eV and E± = 0.61±0.07 eV. The values are 
consistent (within the errors) with the optimum combina- 
tion of energies obtained from fitting average chain length 
dependence on coverage from in-vivo measurements - see 

Fig. EH). 

Further we simulated experimentally obtained scaled 
chain length distributions. Fitting was performed for 
the two in-situ measurements which provided the largest 
data sets (0.04 and 0.08 ML coverage) and for the in-vivo 
measurements. As an example of analyzed data we show 
in Fig. [7] simulations of experimental data for a distri- 
bution function. The simulations were performed for the 
energies E\\ = 0.64 eV and E± = 0.67 eV and correspond 
to a layer with coverage 0.04 ML (see Table HI| . The 
plot, in addition to points representing the whole chain 
population for s > 2, contains the points correspond- 
ing to single indium atoms trapped on C-defects. The 
deviation used for the final fitting was a sum of Dm val- 
ues calculated using the relation Q for the three chosen 
measurements (see Fig. [BJ;). Comparison of experimen- 
tal and simulated data resulted also in almost isotropic 
diffusion of In on Si(100) with activation energies very 
close to the values obtained by the fitting represented by 
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FIG. 7: Experimental and simulated scaled chain length dis- 
tribution functions for a relaxed layer at coverage 0.04 ML 
(in-situ experiment). Results of various simulation experi- 
ments for the same parameters En and E± are plotted. Data 
representing single indium atoms trapped on C-defects are 
plotted as well (on the very left side). 



the diagrams a), b) and d). 

In addition, we investigated how the direction of de- 
tachment of atoms from chains affects values of the es- 
timated diffusion barriers. We found that the direction 
of detachment does not affect results significantly when 
diffusion is nearly isotropic, it plays an important role 
only in case of very anisotropic diffusion. Table HTl shows 
the energy barriers obtained from simulations of different 
experiments, fitted using both the parallel and perpen- 
dicular direction of detachment. Any combination of the 
energies outside the range given by included errors re- 
sults in a double deviation (compared to the best fit) be- 
tween experimental and simulated data. The activation 
energies obtained for the parallel detachment are slightly 
lower than for perpendicular one. 

The diffusion barriers for indium determined from 
KMC simulations (see Table |TT| correspond to almost 
isotropic diffusion. It is in contrast to Albao's et al. 
anisotropic results (En = 0.40 eV; E_l = 0.81 eV). The 
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TABLE II: Activation energies for In adatom diffusion par- 
allel and perpendicular to In chains for two directions of de- 
tachment. Values were obtained by comparing data from two 
in-situ and one in-vivo experiments and simulations for two 
models considering different atom detachment from chains. 

Parallel detachment Perpendicular detachment 
E,| [eV] E x [eV] Ejj [eV] E x [eV] 
0.04 ML 0.64 ± 0.03 0.62 ± 0.07 0.64 ± 0.03 0.67 ± 0.07 
0.08 ML 0.62 ± 0.03 0.61 ± 0.07 0.61 ± 0.03 0.64 ± 0.07 
in-vivo 0.60 ± 0.10 0.65 ± 0.05 0.65 ± 0.10 0.65 ± 0.05 



vales obtained theoretically for Al by Brock and Kelly 8 
are anisotropic as well (Em =0.1 eV; Ej_ = 0.3 eV) but 
much lower. 

Scaled chain length distribution function. It follows 
from the nucleation theor y 20 i 24 that during the submono- 
layer growth the density (per site) N s of islands composed 
of s atoms fulfils the scaling form: 



n s « e 



(i) 



where is coverage and / (s) represents the mean 
island size density. The function f(x), x = s/(s) is the 
scaled island size distribution function. The relation ([T]) 
was confirmed by simulations using different models of ir- 
reversible 2D aggregation and from STM experiments. In 
most cases the shape of the scaled distribution function 
is monomodal - with a peak for x = 1. A monomodal 
function was observed both for homogeneous and het- 
erogeneous nucleation. The scaling for the passage from 
irreversible to reversible aggregation was examined the- 
oretically in Ref. [U The unconventional shape of the 
scaled distribution function - monotonously decreasing 
- observed for ID growth of Ga on Si(100) by Albao 
et al. (and explained by means of KMC simulations using 
strongly anisotropic surface diffusion)- was theoretically 
investigated by Tokar and DreysseiA. They showed that 
for equilibrium growth and a model with atomic inter- 
actions restricted to only nearest neighbors, the scaled 
distribution function is exponentially decreasing. Here 
we obtained for ID submonolayer growth of In similar 
monotonously decreasing chain length distribution func- 
tions - Fig. |1 and Fig. [7J 

If we consider homogeneous nucleation with the de- 
tachment only and the C-defects are excluded from 
our model the simulated distribution function remains 
monotonously decreasing. On the other side excluding 
the process of detachment (irreversible model) results in 
a conventional monomodal distribution function indepen- 
dently on presence of the C-defects in the model. For the 
irreversible model the monotonous distribution can be 
simulated only when strong diffusion anisotropy is intro- 
duced - similarly as reported Albao et al£ 

We conclude that the monotonous form of the chain 
length distribution function obtained for indium layers 
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FIG. 8: Post-deposition relaxation of a simulated scaled dis- 
tribution functions for two detachment barriers - (a) 0.82 eV 
and (b) 0.76 eV. The transition to the monotically decreas- 
ing function corresponding to an equilibrium state is visible. 
The simulation corresponds to a layer with coverage 0.08 ML 
deposited at rate 0.002 ML/s. 



with low coverage (< 0.15 ML) at RT and used depo- 
sition rates can be explained by the process of atom 
detachment from indium chains during the growth. 

Postdeposition relaxation. Further we used our re- 
versible growth model for exploring postdeposition re- 
laxation indicated by experimental data obtained from 
STM measurements at various time after deposition (see 
Fig. [4]). We simulated the growth using the diffusion 
barriers determined and deposition rate 0.002 ML/s for 
two different values of the detachment barrier. The time 
evolution of the distribution functions obtained for vari- 
ous time after deposition is in Fig. [8^, b. If the detach- 
ment rate is small enough with respect to deposition rate 
(a high energy barrier for detachment) the distribution 
function is monomodal just after the deposition and re- 
laxes into a monotonous one - Fig. [8k,. If the detachment 
rate increases the distribution function changes from the 
monomodal to monotonous form even during the growth 
- see Fig. [H]b. The simulation shows that the observed 
system reaches an equilibrium state after ~ 6 hours (es- 
timation for the used experimentally determined barrier 
Edet — 0.82 eV), though the most dramatic change oc- 
curs during the first 10 minutes after the deposition. 

Influence of C-defects. The time constant for the de- 
tachment of an In atom from an adsorption site at a C- 
defect is « 100 x bigger than for the detachment from In 
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FIG. 9: a) Dependence of average chain length on coverage 
simulated for various concentrations of C-defects. Triangles 
represent data from the corresponding in-vivo experiment, 
b) Average chain length is controlled by concentration of C- 
defects. Note also the relaxation of average chain length dur- 
ing time. E\\ = 0.62 eV and E ± = 0.58 eV, E det = 0.82 
eV. 



chains, making the C-defects practically perfect diffusion 
trapsiS. To describe a role of C-defects at heterogeneous 
nucleation and submonolayer growth of In quantitatively 
we simulated growth using the same parameters as for 
a real experiment, only concentration of the C-defects 
was changed. Fig. [9K shows dependence of average chain 
length on coverage obtained for parameters of the in vivo 
experiment with a low deposition rate 1 x 10~ 4 ML/s. 
The average chain length is controlled by concentration 
of the C-defects. The dependence of the average chain 
length on the C-defect concentration is on Fig. It 
can be seen that an effect of postdeposition relaxation 
disappears with the increasing concentration. 



sufficiently high deposition flux the growth of chains 
due to attachment of deposited atoms dominates and 
the effect of the detachment is suppressed. The growth 
becomes irreversible and the scaled distribution function 
has a monomodal form (at almost isotropic surface 
diffusion) . 

Comparison with an equilibrium model. The model of 
Tokar and Drysse^i for equilibrium homogeneous growth 
can provide the only parameter, the nearest-neighbor 
interaction energy V§ N . It can be determined from the 
dependence of average chain length on coverage (using 
the equation (14) in Ref. [TTh . In case of the equilibrium 
model growth characteristics are independent of a 
kinetic path. Our data from in-vivo experiments at a 
rather low deposition rate may reflect a situation not 
too far from an equilibrium state, but the considered 
growth is heterogeneous. The data in Fig. \Qp can be 
approximated by an exponential scaling function. The 
value V§ N — — 0.17 eV obtained from our experimental 
data for indium is similar to the pair-interaction energy 
determined in Ref. [ill for gallium (—0.192 eV). 



Simulations of experimental data showed that the de- 
position rate is a crucial parameter for growing indium 
on Si(100) surface at RT and determines a transition be- 
tween irreversible and reversible character of the growth. 
The measured growth characteristics depend on a process 
of postdeposition relaxation, which has to be included 
into the growth simulation. The monotonous character 
of the scaled chain length distribution function can be 
related to a mechanism of atom detachment from the 
chains. The deposition rate and substrate temperature 
can be used for controlling competition between kinetics 
and equilibrium. The model formulated for simulations 
is too simplified to be used for explanation all experi- 
mentally observed details^ - for example a plateau in 
average chain length dependence on coverage within 0.05 
and 0.12 ML (see Fig. [5]) obtained from two various in- 
vivo measurements - but the model explains the most 
important features of the growth of indium on the sur- 
face Si(100)-2xl at room temperature. 



Growth scenario. The simulations and analysis of 
results allow to suggest a scenario for the experimen- 
tally observed submonolayer growth and discuss a role 
of participating processes at RT. Simulations of the 
postdeposition relaxation show that only a population 
of atoms detached from chains exist on the surface after 
the deposition. They nucleate as new chains or attach to 
other chains and the system moves towards dynamical 
equilibrium. During the deposition, if the deposition 
flux is small enough with respect to detachment rate 
and surface mobility, the growth proceeds at thermal 
equilibrium and a postdeposition effect consists from 
thermal fluctuations only - the scaled chain length 
distribution function is monotonously decreasing. At 



VI. CONCLUSIONS 

STM technique was used for studying growth of indium 
on the Si(100)-2xl surface at low coverage and room 
temperature. Direct observation during the deposition - 
in-vivo measurements - showed that the C-defects act as 
nucleation centers for indium adatoms. The majority of 
indium chain is pinned on one or both ends by a C-defect, 
that determines the average chain length for a given cov- 
erage. The in-vivo observations further revealed the re- 
versible character of the growth due to atom detachment 
from the chains. Statistical characteristics of the In layers 
(average chain length, average length of chains not termi- 
nated at C-defects, dependence of average chain length 
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on coverage, scaled island-size distribution function) were 
obtained from the experiments of two types - the in-vivo 
measurements and the standard in-situ observations af- 
ter the deposition of various coverages. 

The atomistic model with anisotropic diffusion which 
includes presence of C defects on the surface as well as 
detachment of atoms from the chain was developed. Both 
in-vivo and in-situ experiments were simulated using 
KMC method. The simulations showed that the pro- 
cess of atom detachment can explain the monotonously 
decreasing shape of the scaled chain length distribution 
function. Free parameters of the model - activation en- 
ergies for anisotropic surface diffusion - were determined 
by comparison of experimental and simulated character- 
istic of the indium layers. The values obtained for the 
activation energies (see Table [TTJ) correspond to almost 
isotropic surface diffusion in contrary with anisotropic 
data reported for the same group metals Ga and Al ear- 
lier. 
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Appendix 

The deviations between experimental and simulated 
characteristics were determined as follows: 

(i) for the average chain length (s): 

D<i simulated experiment ) ; (2) 

(ii) for the dependence of average length on coverage: 

8=0.25 ML 

D H = V -L(( s )e m " te * ed -(s)se Pe " me " t ) 2 . ( 3 ) 

6=0.025 ML CTe 

where Cg, is a weight parameter - mean square deviation 
of chain length at a given coverage obtained from 
simulation experiments; 

(iii) for the scaled chain length distribution functions: 

S—QO 

jj ^ y ^simulated ^experimental \£) 

s=l 

where function values f s are calculated for each chain 
length of s atoms contained in the data. 
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